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Abstract 

It is proven that the lowest excitations Ei ow {k) of one-dimensional half- 
integer spin generalized Heisenberg models and half-filled extended Hubbard 
models are 7r-periodic functions. For Hubbard models at fractional fillings 
Eiow(k + 2fc/) = Ei ow (k), where 2kf = trn, and n is the number of electrons 
per unit cell. Moreover, if one of the ground states of the system is magnetic 
in the thermodynamic limit, then Ei ow {k) = for any k, so the spectrum 
is gapless at any wave vector. The last statement is true for any integer or 
half-integer value of the spin. 
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Low-energy excitation spectra of correlated systems have long been in the center of the 
theoretical condensed matter research (see |1J and references therein). In 1932 Hans Bethe 
proposed his Ansatz method [[| for correlated wave functions, that was then used to solve 
many one-dimensional quantum and also two-dimensional classical models (see 0], [|J, ||). 
Using the Bethe Ansatz (BA) method des Cloizeaux and Pearson [[]] derived the energy of 
low lying S = 1 (triplet) states of the spin-1/2 antiferromagnetic Heisenberg model (AFHM). 
The spectrum of the excitations was found to be E t (k) = 1/Uir | sin A; |. E t (k) is a tt- 
periodic function of k. des Cloizeaux and Pearson assumed the excitations E t (k) to be the 
elementary excitations of the system, therefore assigning to an elementary excitation spin 1. 
Later investigations |7[ have shown that the spectrum of || was incomplete. The full spec- 
trum of the S = 1 magnon excitations is characterized by two continuous quantum numbers. 
Fixing the wave vector k one still has a continuous one-parametric family of excitations. The 
des Cloizeaux and Pearson spectrum is only a lower bound of these excitations. As a result 
S = 1 states can not be considered as elementary indivisible particles having only one wave 
number k. They are rather two-particle compositions. The problem was finally clarified by 
Faddeev and Takhtajan in Ref. ||. It was shown that all excitations of the spin-1/2 AFHM 
are superpositions of spin-1/2 elementary excitations (called kinks or spinons). For periodic 
rings having odd/even number of sites only odd/even number of kinks are possible in the 
system. Kinks do not form bound states. Their interaction manifests itself only through a 
scattering amplitude. Energies and momenta of many-kink states add, as for independent 
particles. The dispersion relation of a kink is l/2J7rsinfc, with < k < n. The Brillouin 
zone is therefore only half of the original one. The S = 1 states of the model are pairs of 
kinks. The lowest excitations are found by minimizing E (k\) + E (A^) for ki + k 2 = k. The 
result is ki = k and k 2 = (or vice versa), so the lowest excitations probe the one-spinon 
excitations spectrum. For the ferromagnetic case the one magnon spectrum of the one- 
dimensional spin-1/2 ferromagnetic Heisenberg model (FHM) has the form 2J(1 — cos/c), 
possessing a gap at nonzero k and not being 7r-periodic. Careful analysis shows, that one- 
magnon states are not the lowest energy states of the system. Magnons can form bound 
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states and the interaction is so strong,, that the two-magnon states have energy less than 
the energy of one magnon. It was shown, using Bethe Ansatz [[J , that the elementary ex- 
citations of the system (excitations having one quantum number) are magnons and strings. 
Strings are complexes of 2M + 1 spins, having a dispersion relation J 2 m + i (1 ~~ cos A;). The 
lowest bound for the spectrum of excitations is formed, when M goes to infinity. There- 
fore the spectrum is gapless for all wavevectors k. As will be shown below, this property 
characterizes all Heisenberg-like one-dimensional models (integer-spin or half-integer-spin) , 
which are magnetic in the thermodynamic limit. The excitations spectrum of the half-filled 
one-dimensional Hubbard model was derived in [|10f| . The lowest-lying triplet states of the 



model are given by equation (16) of [10]. The function E t (k) is ^-periodic and reaches its 
maximum at q = it/2. For a quarter filling the spectrum of the lowest energy excitations is 
periodic with a period of n/2 = 2kf, where kf is the Fermi wave vector of the noninteracting 



system. For the attractive Hubbard model the excitations spectrum is calculated in [ I I 



The elementary excitations of the model at half-filling are bound pairs and "free" electrons 
(instead of spinons and holons in the repulsive case). The energy- momentum dispersion re- 
lation of excitations is the same, as in the repulsive case, the lowest excitations being again 
7r-periodic. 

It is proven below that the 7r-periodicity of the lowest-lying excitations is a model in- 
dependent feature and holds for all one-dimensional isotropic Heisenberg models having a 
half-integer value of the spin per unit cell, and for all one-dimensional isotropic Hubbard 
models with an odd number of electrons per unit cell. For Hubbard models at fractional 
fillings the spectrum of lowest excitations is periodic with a period im where n is the num- 
ber of electrons per unit cell. For systems with a nonzero groundstate magnetization we are 
able to prove a stronger statement. We call a 1-D system magnetic in the thermodynamic 
limit, if there exist such N , that, for all iV > iV , one has Sn/N > e > 0, where Sn is the 
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spin of the groundstate of a iV-site ring.Q We prove, that if a generalized Heisenberg or 
Hubbard model is magnetic in the thermodynamic limit, then the spectrum of excitations 
is gapless at any wave vector k, i.e., that, for any k, there exists an excitation of arbitrarily 
small energy. The only restriction on the interaction is, that the interaction is short-ranged 
enough, namely, that it falls off more rapidly as 1/r. In 1933 Felix Bloch [12 formulated his 



theorem stating the absence of a stationary current in the ground state of a solid with no 
external field. The physical justification of it is straightforward. Suppose the ground state 
of a solid has got a nonzero current. Let us make a closed electric circuit, consisting of the 
solid and some dissipative device (resistance). Then after a relaxation time the current will 
disappear and some amount of energy will be absorbed in the resistance. Since the energy 
is conserved, the final zero-current state of the solid will have energy less than the initial 



state. A mathematical proof of the Bloch theorem was given by Bohm in ||13fl . The Bloch 
theorem can be also formulated for the spin current, i.e., the current of the ^-component of 
the spin. Both of the versions of the Bloch theorem have been widely used in solid state 
physics, for instance in the theory of superconductivity jnj] . It is proven in this paper, that 
the spin current (current of the z-component of the spin) is zero for Heisenberg and Hubbard 
models in any dimensions. To prove our assertions we use the transformation introduced in 



15fl to prove the vanishing of the excitation gap in the thermodynamic limit of half-integer 



spin Heisenberg models. This transformation is an analogue of the transformation |13| for 
the case of lattice spin models. Consider first a one- dimensional half-integer spin Heisenberg 
model on a periodic ring, having 2N sites. 

H = J ^2 S n • S n+ i — J S*S^ +1 + - ( y S^S n+1 + S^ +1 S n ) (1) 

The sign of J is not specified, therefore both ferromagnetic and antiferromagnetic cases 
are considered. The system is translationary invariant and the wave vector A; is a good 



1 If the groundstate is nontrivially degenerate, then this condition must be satisfied for at least 
one of groundstates. 
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quantum number. For a finite ring, k takes values |^/, where / is an integer. The lowest 
energy excitations of the system are defined as excitations minimizing the energy for a fixed 
value of k. Below we prove that the spectrum of the lowest energy excitations Ei ow (k) is a 
7r-periodic function of k, Ei ow (k) = Ei ow (k + ir). The system has full rotational symmetry, 
therefore it is enough to consider the excitations having — 0- All other excitations can 
be produced from = by applying operators Sf^, ^[ u \\ ■ Consider an eigenstate | ip k ) 
of H, having wave vector k and energy E. 

T\ij k )=e lk \ij k ), H\ij k )=E\ij k ) (2) 

Here T is the translation operator, n — > n + 1. Let us write the state | ip k ) as a linear 
combination of S| u n = spin configurations ip k = J2a A a \ a) . 

Any S^ n Y[ = configuration is characterized by 2NS numbers xi < x 2 < x 3 ... showing 
the positions of fictitious particles, each of them increasing S z on a given site by one. The 
vacuum of the system is assumed to have all spins down. Since the system is periodic, the 
x's are defined only mod 2N. Consider the state 

^ + . = £ew(- + - + ---)A CT |a>. (3) 

cr 

Adding 2N to any of the x's does not change the ip k + n - Therefore the mentioned freedom 
of defining the x's is satisfied. If the operator T acts on the state vp k+w , then all x's are 
incremented by 1. Therefore an additional phase factor e lw2S = e t7T is acquired (2S is odd), 
and the state ip k+7T has a wave vector k + n. In the same fashion the state 

^ k . n = J2e-^ + ^A a \a} (4) 
<j 

has a wave vector k — it. Note that k + 7r and k — ir are the same wave vectors since they 
differ by 2n. Let us evaluate the expectation value (ip k+7r \ H \ ip k+7T ). The Ising part of the 
Heisenberg Hamiltonian does not flip spins. In contrast the exchange part may change the 
position of one spin-up fictitious particle, which amounts to a phase factor e ± ^, depending 
on a direction of the move. This explains the result of a formal calculation, which gives: 



{ljj k+7r | H | Ipk+ir) = (5) 

= E J [1/2 (e*%* I S+S~ +1 | V»*> + e~®(fa | | ^)) + {fa I | . 

To understand the physical meaning of the operator S£S~ +1 — S^ +1 S^ let us find the time 
derivative of the S*. Due to basic principles of quantum mechanics one has |TjJ 

where J* = § (S^S' - S+S'^). Therefore the operator J* = g {S^S' - S+S~_^ 
has a meaning of the current of the z component of the spin through the bond n — l,n. 
The time derivative of is then a difference of incoming and outcoming currents. For any 
eigenstate of the Hamiltonian H all expectation values are stationary, therefore 4+{Sfy = 
{Jn) ~ (J'n+i) = 0- As a result for any eigenstate of H the expectation value {J*) does not 
depend on n. The spin-current operator J n is a pseudovector with components J~£, J*. 
For the states having a definite value of it is 

U?) = (JJ) = 0. (7) 

This is the selection rule for vector operators discussed in §29 of [TH|. Taking a thermody- 
namic limit of (f|) one finds 

{fa +n | H | fa + „) -(^1^1 = 2iV^(e-w(^ | | V*) + e^(^ | S+S~ +1 \ fa) - (8) 

-(V'fc I 5 n+l 5 n I - (^fe I 5 n I V'fc)) = ^ ((^fc | S+ S~ +1 \ fa) - {fa \ S+ +1 S~ \ fa)) . 

In exactly the same way: 

(fa-* | H | fa_^) ~(fa\H \ fa) = -Mi: ({fa | S+S~ +1 | Vfc) - {fa I | fa)) . (9) 

Note that the expressions @, fl9|) differ by sign. Therefore one of the states (Q), (f|) has 
an expectation value of the energy lower or equal than the state | fa). In other words, 
having some excitation at the wave vector k one can always construct an excitation at 
the wave vector k + tc having less or equal energy. Let us now take ipk to be a minimum 
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energy excitation for a given value of the wave vector k. As has been proven there exists 
an excitation with a wave vector k + 7r having the energy less or equal to the energy of fa. 
Therefore: E\ ow (k) > E\ ow [k + 7r). But we could have also started from the state k + 7r to 
show that Ei ow {k + n) > Ei ow (k). Therefore Ei ow (k) = Ei ow {k + n), so the spectrum of 
the lowest energy excitations is 7r-periodic. Equal sign in the inequalities above are realized 
if and only if (J*) = for any state ipi ow (k). Combining this fact with (0) one has J n = 0, 
for the lowest excitations of H having «S| u u = 0- Now consider integer spin models. Then 
the transformations (|f), (f|) produce states with the same wave vector as the initial state. 
Therefore the 7r-periodicity does not work. The statement, that the spin current is equal to 
zero for the lowest-lying = states is still valid (otherwise we would find a state with 
the same k having lower energy). For the spin-1/2 AFHM case it is easy to understand, 
why the spin current is equal to zero for the lowest-lying states in the thermodynamic limit. 
As has been already said, the lowest-lying excitations are two-spinon excitations, so they 
involve two spin-1/2 particles. To produce a nonzero spin current through each bond in 
the thermodynamic limit one would need to excite a finite number of the particles per unit 
volume(length). 

Inclusion of more than nearest-neighbor interactions or next order terms like (S„ • S„ + i) 2 
does not change anything in the proof. Let us, for instance, add to the Hamiltonian ([!]) a 
next-nearest-neighbor interaction term a (S n ■ S n+ 2). Here a is a coupling constant. Then 
the equations (@)-([|) are not changed. The right-hand side of the equation (|5[) acquires an 
additional term 

a £ [1/2 (e 2 w<^ | S+S' +2 \ Vfc) + e~ 2 w(^ | S+ +2 S~ \ fa)) + (fa \ S^ +2 S z n \ fa)} . (10) 

n 

This term is obtained by a direct calculation, and its physical explanation is similar 
to the explanation of the equation The Ising part of the next-nearest-neighbor term 
a (S n • S n+ 2) does not flip spins. In contrast to that, the exchange part may change the 
position of one spin-up fictitious particle by two lattice sites, which amounts to a phase 
e ±2 w, depending on a direction of the move. Taking the thermodynamic limit, one has an 



additional term 

2ia-K ((V> fe | S+S~ +2 | ip k ) - (ip k | S+ +2 S~ | xpk}) . (11) 

on the right-hand side of the equation @. In the same fashion the equation (Q) acquires 
the term 

- 2iair ((ip k | S+S~ +2 | ip k ) - (ip k \ S+ +2 S~ | V*)) • (12) 

Then arguments, identical to those given after equation @, lead to a 7r-periodicity of lowest- 
lying excitations. One can now proceed adding next-next-nearest-neighbor interactions, etc. 
The proof goes exactly in the same way, as soon as the interaction range is finite. If the 
interaction has a long range behavior, then one must require the long-range part to fall off 
quickly enough. Namely, assuming that the interaction falls off as l/r 1+7 , 7 > 0, one can 
successfully repeat the steps to prove the 7r-periodicity conjecture in this case. I 

leave a detailed proof for this case to a longer paper ||17|| . One can also add nonlinear terms 
like (S n • S n+ i) m , and, repeating the same steps, show by a direct calculation, that the 7r- 
periodicity conjecture holds for this case too. Let us, for instance, add to the Hamiltonian 
([!]) a term (3 (S n • S n+ i) 2 . We can rewrite this term as 
This expression contains terms which do not change Xj's, terms that change one of x^'s by 
1, and terms, like (5'^ +1 S'~) 2 , that change positions of two Xj's by one. An interested reader 
might check again, that the steps (0)-© can be straightforwardly repeated. 

To formulate it shortly, for any one-dimensional half-integer-spin per unit cell isotropic 
Heisenberg Hamiltonian the 7r-periodicity conjecture is true. The proof of the general con- 
jecture goes exactly in the same way, as the examples considered above, and will be given in 
the longer paper ||17|| . One can also note, that the full spherical symmetry is not necessary. 
The only facts really used are, that the Hamiltonian is symmetric with respect to a rotation 
around the z-axis, and that the lowest energy excitations can be chosen to have = 0- 
As long as these conditions are satisfied, the spectrum of the lowest energy excitations is 
7r-periodic. 
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As have already been stated, the proof works also for long-range interactions, if the in- 
teraction goes to zero more rapidly than 1/r. Haldane and Shastry introduced an integrable 
Heisenberg model with an interaction going as 1/r 2 . The spectrum of lowest spinon excita- 



tions of this model is given in and is 7r-periodic, which is an example of the 7r-periodicity 
conjecture, considered here. 

Let us note, that in some models not only the lowest excitations might be 7r-periodic, but 
the groundstate itself can break the discrete translational symmetry and dimerize. In this 
case the size of the unit cell is doubled, and the Brillouin zone is halved. Then all the exci- 
tations of the system will possess a 7r-periodicity, which is in this case related to a symmetry 
breaking and dimerization. An example of such a system is a frustrated 1-D spin-one-half 
Heisenberg chain with J 2 ^> J\ , where J\, J2 are antiferromagnetic nearest-neighbor and 



next-nearest-neighbor couplings. Extensive studies of this model were performed in p9 
This model may also be viewed as a zigzag chain, with an intrachain coupling J2 and an 
interchain coupling J x - The model is exactly solvable for J 2 = 1/2 Ji [pm , the groundstate 
is dimerized, and the discrete translational symmetry is broken. Density Matrix Renormal- 
ization Group studies ( []TU[) show, that the discrete symmetry breaking appears not only in 
the exactly solvable case, but for all J2 > 0.25Ji. There are also helical correlations in the 
model, which could lead to an incommensurate order in higher- dimensions, but which are 
destroyed in one dimension, leading to a finite correlation length and to a gapped state. 

Let us now turn to Hubbard models. In this case we claim the following 2/cj-periodicity 
conjecture to be true. For any isotropic one- dimensional extended Hubbard model the 
spectrum of lowest-lying excitations is 2/c/-periodic, where 2kf = nn mod 2tt, and n is the 
number of electrons per unit cell. Let us note, that for many simple one-dimensional systems 
this definition of kf coincides with the Fermi wave vector, when the hubbard U is switched 
off, but for more complicated systems (for instance for systems with partially filled bands), 
it is not the case. Consider first the simplest case of a one-band Hubbard model 

H = t £(4^+1.- + h.c.) + U £(n i|T n u ) • (13) 



having n electrons per site, n is a rational number, n = C1/C2, where c\ and c 2 are integers. 
We consider a ring of IciN sites, eventually taking a limit iV — > oo. The number of electrons 
on the ring is equal to 2c\N and is an even number, therefore the full spin of the system 
takes integer values. Because of the full rotational symmetry, it is enough to consider only 
the states ^\ = 0- All other excitations can be produced from = by applying 

operators Sj-^jj' • Analogously to @ we consider an eigenstate | ipk) of H, having a 
wave vector k and energy E. We write the state | ipk) as a linear combination of = 
electron configurations ipk = J2a A a \ a). Each of the configurations | a) is characterized by 
ciN numbers x\ < X2 < 23... , showing the positions of C\N spin-up electrons, and by c\N 
numbers x[ < x' 2 < x f 3 ..., showing the positions of spin-down electrons. Since the system is 
periodic, these positions are defined only mod 2Nc2- Now, in analogy to @, consider the 
state ipk+2k f = Ha e 2JVc2 ( X1+X2+ -} | a y Adding 2Nc2 to any of the x's does not change 
the ipk+-K- Therefore the mentioned freedom of defining the x's is satisfied. If the operator 
T acts on the state ipk+ir, then all x's are incremented by 1. Therefore an additional phase 
factor e i7rCl//c2 = e t2k f is acquired, and the state ipk+2k f has a wave vector k + 2k f. Now 
we need to construct a state, analogous to the state (^). Since k + 2k f and k — 2k f are 
in general different wave vectors, simply taking the state k — 2 kf will not work. But we 
note, that since n is rational, there always exist an integer number X, such that k — X2kf 
and k + 2kf are the same wave vectors. Let us take X = 2c2 — 1. Then k — X2kf = 
k — (2c2 — l)7rc 1 /c 2 = k — 2ttci + txci/c2 = k + 7rci/c2 = k + 2k f. Consider now the state 
ipk-X2k f = Ha e X 2Nc2 ^ 1+X2+ '"- ) A a I a). This state has a wave vector k — X2kf, which is, as 
we already said, the same as k + 2k f. 

Let us now evaluate the expectation value (ipk+2k f \ H \ ipk+2k f )- The interaction U, and 
the spin-down hopping do not change the positions of spin-up electrons. In contrast the 
spin-up hopping may change the position of one spin-up electron by one, which amounts to 

I 2tyi 

a phase factor e 2c ^ N , depending on a direction of the move. This explains the result of a 
direct calculation, which gives, that in the thermodynamic limit iV — > 00 one has 
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(*Pk+2k f I H I ipk+2k f ) ~ (i>k I H I ipk) == it-it (j^ k | c+ T c n+ i iT | ^>fc) - I c n+l,T C n,T I Vte)) • ( 14 ) 

{i>k-X2k s I # I ipk-X2k f ) - {^k\H\ ip k ) == -Xit-K (jipk I c^ T c n+ i >T | Vfc) - (ipk I c+ +liT c„ iT I tpk}) ■ 

We see, that, since X is by construction a positive integer, if one of the above expressions 
is positive, the other is negative, and vice versa. Therefore one of the above states has the 
expectation value of energy less or equal to the energy of the state | ipk)- Note that both of 
considered states have the wave vector k + 2kf. Now repeating arguments given after the 
expression @, we state that the lowest excitations are periodic functions with a period 2k f. 
The 2fc/-conjecture is therefore proven. It is again possible to check, that adding more than 
nearest-neighbor hoppings does not change the picture. One can also add arbitrary intersite 
interaction terms, of the type VkUin^k, since these terms are diagonal in the considered 
basis and drop out of calculations. One can also consider 1-D lattices with more than one 
site per unit cell (for instance n-leg ladders), to show again, that the spectrum of excitations 
is 2k f periodic, where 2kf we define as 7m, n being the number of electrons per unit cell. 

So far we used for our derivations the subspace — 0. Consider now a Heisenberg 
model, where one of the (possibly many) ground states is magnetic in the thermodynamic 
limit, i.e., limTv^oo 5 , f u u/2iV > M > 0, where S^}} is the spin of the ground state. Then 
for the fixed iV the °f the ground state can be chosen to be from —2MN to 2MN. 
Let us choose the ground state to have Sf^v] = C, and write it in a linear combination of 
'"'full = ^ con fig ura ti° ns - Let us now consider the state | ip') = J2a e G ^( Xl+X2+ "^A a | a) , 
where xi, X2, ... are again the positions of the + 2C fictitious particles. G is an integer 
number. The wave vector of the state > differs from the wave vector of the ground 
state by 2ttG(S + C/N). Now choosing different (allowed) values of C and G, one can show 
after a careful analysis (which will be given in details in [|I7]]), that in the thermodynamic 
limit there is no gap at any value of k. This result does not depend on whether the spin is 
half-integer or integer. Therefore if the system is magnetic in the thermodynamic limit, the 
spectrum is gapless for any k. A similar result can be obtained for one-dimensional extended 
Hubbard models Hl7f . 
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